MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS  1963  A 


v  v  vlv  *  •**  ■*  *  -*• 


’  \*  V  •  •  V  V  V 

■  •  v  ‘A-  v  vVj-  v 

“  -  ■  -  'j-  • 


■  .'•.'•■.'•'A".'''. v'  .".-  '  .’  - .  \.'  •  •  .  • 


AD- A 171  764 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data^Entered)*  "  ’ 

REPORT  DOCUMENTATION  PAGE 

'  READ  INSTRUCTIONS  .S 

BEFORE  COMPLETING  FORM 

i.  report  number 

AFIT/CI/NR  86-151T 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT'S  CATALOG  NUMBER 

4  TITLE  (and  Subtitle) 

Improving  Risk  Characterizations  Based 

On  Time  To  Response 

S.  TYPE  OF  REPORT  &  PERIOD  COVERED 

THESIS 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AuTHORfx) 

Ronald  J.  Berdine 

8.  CONTRACT  OR  GRANT  NUMBE Rtf's) 

9  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

AFIT  STUDENT  AT:  Texas  A&M  University 

to.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 

CONTROLLING  OFFICE  NAME  AND  ADDRESS 

12.  report  DATE 

1986. 

13.  NUMBER  OF  PAGES 

69 

MONITORING  AGENCY  NAME  ft  ADDRESSfll  different  from  Controlling  Olfice) 

15.  SECURITY  CLASS,  (of  this  report) 

UNCLASS 

15a.  DECL  ASSI  F|  CATION  DOWNGRADING 
SCHEDULE 

distribution  statement  (of  this  Report)  |  y  |  '  |  •  | 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED  ELECTE 

DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  If  different  from  Report) 


SEP  1  7  1986 


18  SUPPLEMENTARY  notes 

APPROVED  FOR  PUBLIC  RELEASE:  IAW  AFR  190-1 


YW  E.  WOLAVER 


L'Ttift  E.  WOCAVER  '  StyJfc 
Dean  for  Research  and  9 
Professional  Development 
AF1T/NR 


*9  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  Identify  by  block  number) 


I  20  ABSTRACT  (Continue  on  reverse  side  If  necessary  and  Identify  hy  block  number) 


ATTACHED 


/ ( 


ABSTRACT 


Improving  Risk  Characterizations  Based  on  Time  to  Response 

(Hay  1986) 

Ronald  James  Berdine,  B.  S.,  Iowa  State  University; 

M.  S.,  Stanford  University 

Chairman  of  Advisory  Committee:  Dr.  Robert  L.  Sielken  Jr. 


One  of  the  important  aspects  of  quantitative  cancer  risk  assessment 

is  to  model  the  frequency  of  a  carcinogenic  response  as  a  function  of 

exposure.  Most  of  the  current  time-to-response  modeling  has  been  done 
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using  the  ’’simple"'  model  which  does  not  include  the  effect  of  competing 
risks.  As  a  result,  the  probabilities  which  are  stated  in  terms  of  these 
models  do  not  correspond  to  real  world  occurrences.  Rather  they 
correspond  to  a  fictional  world  where  no  competing  risks  exist;  that  is, 
a  world  where  only  the  risk  being  modeled  is  present.  The  cause-specific 
family  of  models  is  preferred  since  it  incorporates  the  effect  of 
competing  risks.  A  method  is  presented  for  utilizing  life  table  data  and 
existing  simple  model  maximum  likelihood  estimates  to  obtain 
cause-specific  risk  characterizations. 

Representations  for  various  risk  characterizations  (e.g., 
time-to-response  probability,  mean  response  free  period,  mean  free  dose, 
and  virtually  safe  dose)  under  both  models  are  compared.  The  results 
show  that  the  simple  model  always  overstates  the  effect  of  the 
carcinogen.  The  overstatement  can  easily  be  in  the  20 %  -  200 %  range. 


Two  alternatives  are  Introduced  for  modeling  time  to  death  from 
tumor  when  the  cause  of  death  is  uncertain.  One  alternative  is  to  model 
the  time  to  death  Irrespective  of  cause.  The  other  alternative  focuses 
on  the  increase  in  the  hazard  rate  due  to  the  presence  of  a  tumor.  Thus 
deaths  caused  by  this  increase  in  hazard  rates  are  considered  as  being 
related  to  the  introduction  of  the  carcinogen.  Both  alternatives  give 
risk  characterizations  which  may  be  more  relevant  to  cancer  risk 
assessment . 
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is  to  model  the  frequency  of  a  carcinogenic  response  as  a  function  of 
exposure.  Most  of  the  current  time-to-response  modeling  has  been  done 
using  the  "simple"  model  which  does  not  include  the  effect  of  competing 
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family  of  models  is  preferred  since  it  incorporates  the  effect  of 
competing  risks.  A  method  is  presented  for  utilizing  life  table  data  and 
existing  simple  model  maximum  likelihood  estimates  to  obtain 
cause-specific  risk  characterizations. 

Representations  for  various  risk  characterizations  (e.g., 
time-to-response  probability,  mean  response  free  period,  mean  free  dose, 
and  virtually  safe  dose)  under  both  models  are  compared.  The  results 
show  that  the  simple  model  always  overstates  the  effect  of  the 
carcinogen.  The  overstatement  can  easily  be  in  the  20%  -  200%  range. 


Two  alternatives  are  introduced  for  modeling  time  to  death  from 
tumor  when  the  cause  of  death  is  uncertain.  One  alternative  is  to  model 
the  time  to  death  irrespective  of  cause.  The  other  alternative  focuses 
on  the  increase  in  the  hazard  rate  due  to  the  presence  of  a  tumor.  Thus, 
deaths  caused  by  this  increase  in  hazard  rates  are  considered  as  being 
related  to  the  introduction  of  the  carcinogen.  Both  alternatives  give 
risk  characterizations  which  may  be  more  relevant  to  cancer  risk 
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1.  INTRODUCTION 

One  of  the  important  aspects  of  quantitative  cancer  risk  assessment 
is  to  model  the  frequency  of  a  carcinogenic  response  as  a  function  of 
exposure.  The  type  of  modeling  has  generally  been  of  the  form  where  the 
frequency  of  response  is  modeled  as  a  function  of  dose  alone  as  opposed 
to  including  both  dose  and  time.  Only  recently  have  researchers 
attempted  to  use  a  time-  to-response  variable  to  model  more  completely 
the  biological  development  of  tumors  and  associated  responses  (e.g.. 
Society  of  Toxicology  EDqi  Task  Force  (1981)).  The  response  of  concern 
could  be  death,  death  caused  by  tumor,  death  with  tumor,  tumor  onset  or 
any  number  of  other  possibilities.  This  paper  is  concerned  with  effects 
resulting  in  observable  times  to  response.  Methods  for  modeling 
unobservable  times  are  given  by  Mcknight  and  Crowley  (1984)  and 
Kalbfleisch,  Krewski  and  Van  Ryzin  (1983). 

Heretofore,  risk  characterizations  based  on  time  to  response  have 
been  based  on  the  simple  model  where  the  general  form  of  the  probability 
distribution  for  the  time  to  a  specified  response  is 

Ps(t;d)  =  Pr(T£t;d) 

with  d  being  the  dose  and  T  being  the  time  to  the  response  of  interest. 
However,  P$(t;d)  does  not  reflect  the  fact  that  a  competing  response 
might  occur  before  the  specified  response  of  interest. 

Two  models  which  incorporate  times  to  competing  responses  are  the 
latent  failure  time  model  and  the  cause-specific  model.  In  Section  2 

Citations  will  follow  the  format  of  Biometrics. 
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these  models  are  discussed  and  the  cause-specific  model  is  shown  to  be 
the  model  of  choice.  The  general  form  for  the  cause-specific  model  is 

Pcs(t,j;d)  =  Pr(T<t,J=j;d) 

where  there  are  p  mutually  exclusive  and  collectively  exhaustive 
responses,  J  is  the  indicator  of  the  one  response  observed,  T  is  the  time 
at  which  the  observation  is  made  and  d  is  the  dose. 

In  Section  3  the  formulations  of  both  the  simple  and  cause-specific 
models  and  their  likelihood  functions  are  given.  A  relationship  exists 
between  the  simple  and  cause-specific  likelihood  function  that  allows  new 
risk  characterizations  to  be  calculated  using  results  from  an  experiment 
where  the  simple  model  has  been  applied.  Thus  improved  risk 
characterizations  can  easily  be  calculated. 

Section  4  is  concerned  with  the  selection  of  hazard  functions.  In 
the  cause-specific  model  a  hazard  function  for  each  cause  of  an 
observation  can  be  modeled.  Most  of  the  results  developed  do  not  depend 
on  the  parametric  form  of  the  specified  response's  hazard  function; 
however,  a  very  general  form  is  used  to  generate  examples  and  simulate 
data.  The  competing  risks'  hazard  function  can  be  derived  from  life 
table  data.  For  many  species  such  data  are  available. 

Sections  5  and  6  are  concerned  with  the  impact  on  certain  risk 
characterizations  of  using  the  cause-specific  model  instead  of  the  simple 
model.  The  magnitude  of  this  impact  is  measured  in  terms  of  relative 
bias.  Relative  bias  is  the  fraction  by  which  the  simple  risk 
characterization  overstates  or  understates  a  corresponding  cause-specific 
risk  characterization.  In  certain  cases,  bounds  are  derived  for  the 
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relative  bias  as  well  as  examples  demonstrating  the  actual  relative 
biases.  The  risk  characterizations  to  be  investigated  are  the 
probability  that  the  specified  response  occurs  before  a  time  t,  the  mean 
response  free  period,  the  mean  response  free  dose  and  the  virtually  safe 
dose.  Implementation  of  the  cause-specific  model  requires  that  the 
definitions  of  these  risk  characterizations  be  refined  in  order  to 
incorporate  the  impact  of  competing  risks. 

In  Section  7  two  new  options  for  obtaining  risk  characterizations 
are  investigated.  These  options  are  for  the  special  case  where  the 
original  response  of  interest  is  death  from  tumor  but  the  cause  of  death 
is  uncertain.  Both  of  the  new  options  rely  on  a  new  definition  for  the 
response  of  interest.  One  option  defines  the  specified  response  to  be 
death  from  all  causes  combined  while  the  second  option  defines  the 
specified  -response  to  be  death  related  to  tumor  as  opposed  to  death  from 
tumor. 

Finally,  Section  8  gives  the  conclusions  that  can  be  drawn  from  the 
results  in  Sections  2  through  7.  A  short  discussion  as  to  the  utility  of 
the  procedures  derived  in  this  dissertation  and  their  impact  on  future 
risk  assessment  procedures  is  also  given. 
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2.  REVIEW  OF  THE  LITERATURE 


The  theory  of  competing  risks  was  formalized  by  Chiang  (1970)  and 
Moeschberger  and  Oavid  (1971).  Both  used  a  model  which  is  referred  to  as 
the  latent  failure  time  model.  This  model  uses  p  times  to  failure,  T^ , 
...»  Tp,  each  denoting  the  random  "latent"  time  to  failure  by  one  of  p 
competing  risks.  For  each  individual  only  the  minimum  of  T^,  ...,  Tp 
is  observed.  Thus,  if  T  =  mi n(Tj , . . . ,T  ) ,  then  probability  statements 
giving  the  probability  of  survival  until  a  time  t  are  of  the  form 

Pr(T>t). 

Chiang  (1970)  defines  the  crude  survival  function  as 
0(tjf...,tp)  =  Pr(T1>t1,...,Tp>tp)  , 

and  the  hazard  function  for  each  of  the  latent  times  to  failure  can  be 
defined  by 

Pr(t<T.<t+At | T>t ) 

X  -  ( t )  =  lim  - - - ! - 

J  At-*-0  At 

=  _ll  (log[0( t,  ....  ,t  ) ]}  for  j  =  1,  ...,  p. ( 2 . 1 ) 

dtj  p  t1=...=tp=t 

Also  defined  are  the  marginal  survival  functions 
O.(tj)  -  Pr(T1>0,...,Tj.>tj  ,...,Tp>0)  . 

With  the  assumption  that  T^,  ....  Tp  are  independent  random  variables 
equation  (2.1)  becomes 
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for  j  =  1,  ,  p. 

This  is  the  assumption  used  throughout  the  work  by  Chiang  (1970)  and 

Moeschberger  and  David  (1971).  In  fact,  Tsiatis  (1975)  showed  that  the 

latent  failure  time  model  is  nonidentifiable  unless  the  independence 

assumption  is  made.  Examples  where  more  than  one  model  lead  to  the  same 

set  of  hazard  functions  when  the  independence  assumption  is  not  assumed 

are  found  in  both  Tsiatis  (1975)  and  Kalbfleisch  and  Prentice  (1980). 

Prentice  et  al .  (1978)  discussed  the  latent  failure  time  model 

before  developing  the  cause-specific  model.  They  noted  that  besides  the 

assumption  of  statistical  independence  there  is  another 

"...very  strong  assumption  that  the  time  to  failure  from  cause 
j  under  one  set  of  study  conditions  in  which  all  ...  causes  are 
operative  is  precisely  the  same  as  under  an  altered  set  of 
conditions  in  which  all  causes  except  the  j-th  cause  have  been 
removed." 

Thus,  not  only  is  statistical  independence  required  for  identifiabilty 
but  a  physical  independence  among  causes  of  death  is  implied. 

As  an  alternative,  Prentice  et  al .  (1978)  formulate  the 
cause-specific  model.  It  does  not  require  statistical  independence  nor 
does  it  imply  any  specified  physical  relationship  among  the  causes  of 
death.  (The  cause-specific  model  formulation  is  given  in  Section  3.2.) 
Moreover,  as  the  likelihood  was  developed,  they  noticed  that  it  could  be 
partitioned  into  multiplicative  factors,  each  depending  on  a  set  of 
mutually  exclusive  parameters.  This  property  is  important  in  the 
development  of  methods  used  in  this  dissertation  and  is  discussed  more 
fully  in  Section  3. 
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The  choice  of  a  parametric  form  for  each  cause-specific  hazard 
function  has  received  much  attention  in  the  literature.  Not  only  must 
the  resulting  parametric  model  fit  the  data  well,  but  it  must  provide  an 
accurate  representation  of  the  carcinogenic  effect  at  doses  where  no  data 
may  be  available.  Armitage  (1982)  discusses  the  implications  that 
various  parametric  models  have  on  risk  characterizations  in  the  low  dose 
region.  A  detailed  examination  of  the  effect  of  the  models  on  a 
specified  risk  characterization  (i.e.,  a  virtually  safe  dose  which  is 
discussed  in  Section  6)  is  given  in  Krewski  et  al .  (1983).  One  model 
which  allows  considerable  flexibility  is  given  by  Hartley  and  Sielken 
(1977)  and  Hartley,  Tolley  and  Sielken  (1981).  This  is  the 
Hartl ey-Si el  ken  model  used  in  Section  4  when  a  parametric  form  for  the 
hazard  function  for  the  response  of  interest  is  specified.  A  review  of 
several  parametric  models  is  given  in  Kalbfleisch  et  al .  (1983). 

Some  research  has  also  been  done  on  estimating  tumor  prevalence  when 
the  time  to  tumor  onset  is  unobservable  (e.g.,  Dinse  (1985)  and  Dewanji 
and  Kalbfleisch  (1985)).  Mcknight  and  Crowley  (1984)  showed  that  the 
tumor  incidence  rate  for  unobservable  time  to  tumor  onset  cannot  be 
estimated  unless  interim  sacrifices  are  made.  All  of  the  models  used  in 
the  research  cited  here  use  an  extension  of  the  cause-specific  model  as 
opposed  to  the  latent  failure  time  model. 


3.  TIME-TO-RESPONSE  MODELS 


3.1  Simple  Time-to-Response  Models 

Simple  time-to-response  models  focus  on  the  time  to  a  specified 
response  and  consider  the  times  to  competing  responses  solely  as 
censoring  times  and  not  part  of  the  model.  For  example,  let  the 
specified  response  be  time  to  death  due  to  a  tumor.  If  a  subject  dies 
from  some  cause  other  than  tumor,  then  the  time  of  death  is  treated  as  a 
censoring  time  for  the  time  to  death  from  tumor.  In  the  simple  model  the 
probability  refers  to  the  anticipated  time  of  death  from  tumor  regardless 
of  whether  a  competing  cause  of  death  might  occur  prior  to  that  time  and 
the  death  from  tumor  not  be  observed. 

3.1.1  Model  Formulation 

Suppose  an  experiment  is  conducted  where  several  doses  of  a 
carcinogen  are  administered  to  N  test  subjects.  The  observations  are 
(T.j  ,d.j ) ,  i  =  1,  ....  N  where,  for  the  i-th  subject 

Ti  =  time  to  specified  response  or  censoring  time, 

/* 

1  if  the  specified  response  is  observed  at 

6i  = 

0  if  the  specified  response  is  censored. 
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The  hazard  function  X(t;d)  is  defined  by 

X  ( t ;  d )  =  lim  -P^t<Tl^.At;dlT>t^ 
At-K)  At 


_  Oog[Pr(T>t)l)  , 
dt 


so  that 


P$(t;d)  =  P r ( T_<t ; d )  =  1  -  exp[-Jg  X(u;d)  du] 


=  /g  X(u;d)  exp[-A(u;d)]  du 


(3.1) 


where  A(u;d)  =  /g  X(v;d)  dv  is  called  the  cumulative  hazard  function. 

The  model  given  by  (3.1)  will  be  referred  to  as  the  "simple"  model 
herein.  Each  model  is  uniquely  determined  by  its  hazard  function,  or 
equivalently,  its  cumulative  hazard  function. 

3.1.2  Likelihood 

For  the  i-th  individual,  the  observation  can  be  either  the  response 
of  interest  or  a  censoring  time  such  as  sacrifice  or  death  from  some 
other  cause.  The  possible  contributions  to  the  likelihood  for  the  i-th 
individual  are  given  in  Table  1. 


Table  1. 

Simple  model  likelihood  contributions 
Likelihood  contribution  Observation 


X(ti ;di )exp[-A(ti  ;di )] 

Response  of  interest 

exp[-A(t •  ;d • ) ] 

Censoring  time 
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Let  S  be  the  vector  of  parameters  associated  with  Mt;d).  Also, 
without  loss  of  generality,  assume  that  the  first  observations  are 
the  observed  responses  of  interest  and  the  last  observations  are 
the  censored  observations  where  N  =  ni  +  Then  the  likelihood 
function  is  given  by 

nl  N 

L(©;t)  =  n  { X ( t i  *, d  - )  exp[-A(ti  ;di )]}  n  {exp[-A(t .  ;di ) ] } 
i=l  i=n.+l 

which  can  be  written  as 

N  S>. 

L(6;t)  =  n  {[X(t. ;di )1  exp[-A(ti ;di )]}  . 

7*1 

3.2  Cause-Specific  Models 

In  the  formulation  of  the  simple  models,  the  only  random  variable  is 
the  time  to  the  specified  response.  Thus,  P s ( t ; d )  gives  a  hypothetical 
probability  of  the  specified  response  by  time  t.  The  actual  frequency  of 
the  specified  response  is  never  greater  than  Ps(t;d)  but  can  be  less  than 
Ps(t;d)  due  to  competing  risk  deaths  prior  to  the  time  of  the  specified 
response.  The  cause-specific  model  explicitly  includes  death  due  to 
competing  risks  as  well  as  the  time  to  the  specified  response. 

3.2.1  Model  Formulation 

For  the  i-th  individual,  let  T^  be  the  time  to  one  of  p+1  mutually 
exclusive  responses.  Let  J  =  0,  ...,  p  be  the  indicator  of  the  one 
response  which  is  observed.  Let  J  =  0  correspond  to  a  scheduled 
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sacrifice  and  let  J  =  1  correspond  to  the  specified  response.  Often  it 
is  not  necessary  to  distinguish  among  the  competing  risks,  so  that  p  =  2 
with  J  =  2  corresponding  to  a  death  or  censoring  by  any  competing  risk. 


Also,  let  di  be  the  dose  and  let 

1  if  the  i-th  individual's  response  is  response  j, 
j  =  1.  ....  P 


0  otherwise  . 

Define  the  cause-specific  hazard  function  as 

\ . (t;d)  =  lim 
J  At*0  At 

«  ±  {log[Pr(T>t,J«j;d)]}  . 
dt 

The  cause-specific  model  is  then  given  by 

Prc(t,j;d)  =  Pr(T<t  ,J=j  ;d)  =  /J  X  -  (u  ;d)  exp[-  f  A.(u;d)]  du  . 
cs  “  u  J  k=l  K 

3.2.2  Likelihood 

Kalbfleisch  and  Prentice  (1980)  give  the  development  of  the 
likelihood  for  the  cause-specific  model.  The  contributions  are  as 
follows  (Table  2): 


Table  2. 


Cause-specific  likelihood  contributions 
Likelihood  contribution  Observation 


X j ( t i ; d i )  exp[-  I  Ak(t - ;di )1 
k  1 

exp[-  f  A.  (t.  ;d . )  ] 
k=l  K  1  1 


Response  of  type  j 
j  -  1  *  •••>  P 

Sacrifice  (j  =  0) 


Notice  that  the  i-th  contribution  can  be  written  as 


p  6.  . 

{[^j  (ti  ;d. )]  1J  exp[-Aj.(ti  ;di )]}  _ 


Let  0^,  ....  0^  be  the  vectors  of  parameters  associated  with  A^(t;d), 
...»  Ap(t;d)  respectively.  Then  the  likelihood  function  is 


N  P  6.  . 

L(0.,  0  ;t)  *  H  {H  l>,(t.  ;d. )]  J  exp[- A  ( t .  ;d . ) 1 } 

~P  ~  i=l  j=l  J  1  1  JIT 


P  N  6.  . 

=  n  f  n  [\  (t.;d,)]  J  exp[-A, (t . ;d . ) ] f 
j=l  i=l  J  1  ITT 


P 

=  n  L  •  ( 0  . : t )  (3.2) 

j=l  J  ~J  ~ 


where  Lj (Sj ;t )  is  the  same  likelihood  function  from  a  simple  model  with 
the  j-th  response  being  the  specified  response  of  interest  and  all  other 


responses  are  treated  -s  censors. 
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3.2.3  Properties 

The  cause-specific  model  is  an  intuitively  appealing  model  because 
of  its  ability  to  incorporate  competing  risks.  It  also  has  several 
desirable  properties.  Naturally,  J  and  T  are  not  necesssarily 
independent.  The  simple  model  is  also  a  special  case  of  the  cause- 
specific  model  if  p  =  1.  It  is  very  general  since  the  time  T  does  not 
have  to  be  represented  as  the  minimum  of  p  independent  response  times  as 
in  the  latent  failure  time  model. 

The  tractability  of  the  cause-specific  models  is  facilitated  by  two 
characteristics. 

(i)  In  the  likelihood  function  (3.2),  the  9.,  ...,  6  appear  in 
separate  multiplicative  terms. 

(ii)  Each  of  the  multiplicative  terms  in  the  likelihood  (3.2)  has 
the  same  form  as  the  likelihood  corresponding  to  a  simple 
model . 

The  maximum  likelihood  estimates  of  the  cause-specific  model 
parameters  in  X1(t;d)  are  the  same  estimates  that  would  be  obtained 
using  a  simple  model  and  treating  individuals  who  succumb  to  a  competing 
risk  as  censored.  Therefore,  if  these  simple  model  estimates  are 
available  and  X?(t;d)  is  known,  the  cause-specific  probabilities  and 
other  associated  risk  characterizations  can  be  calculated  without  making 
any  alterations  in  the  estimates  of  the  parameters  in  \j(t;d). 
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4.  CHOICE  OF  HAZARD  FUNCTIONS 

4.1  A  Form  for  the  Hazard  Function  of  the  Specified  Response 

A  parametric  form  for  X^(t;d)  will  be  used  whenever  specific 
examples  or  simulations  are  required.  Theoretically,  any  form  of  X^(t;d) 
could  be  used.  A  special  case  of  the  Hartley-Si  el  ken  model  will  be  used 
with  Aj(t;d)  =  H(t)g(d)  and  the  dose  function,  g(d),  being  a  non¬ 
negative  convex  function  of  dose.  In  particular,  a  Weibull  form  for  H(t) 

1/  p 

will  be  used;  namely  H(t)  =  pt  and  g(d)  =  atg  +  a^d  for  s  >_  1. 

4.2  A  Form  for  the  Hazard  Function  of  the  Competing  Risks 

The  form  for  \2(t;d)  will  be  defined  in  terms  of  S,,(t;d)  = 
expC-AgUjd)].  It  is  assumed  that  the  hazard  function  for  the  competing 
risks  is  not  dependent  on  dose.  Therefore,  the  "d"  notation  will  be 
suppressed  from  this  point  on  when  referring  to  and  S t ) .  Let 

S2(t)  be  a  piecewise  linear,  continuous  function  not  involving  dose  with 

S2(t)  =  a.jt  +  bi  for  ti  £  t  <_  t-+1,  i  =  0,  ....  m  ,  (4.1) 

with  a.j  <_  0,  bi  >  0. 

This  implies 

A2(t)  =  -log  (a.jt  +  bi )  for  ti  <_  t  _<  ti+1  ,  i  =  0,  ....  m  . 

If  the  time  is  scaled  such  that  t.+^  -  t-  =  1,  then  ai  = 

S2 ( t i  +i )  _  S2(ti ).  The  hazard  function,  can  then  be  expressed  as 


Ut)  -  ±  [a  (t)]  ■ 

1  dt  ^  a*t  +  b. 


[s2(ti+i)  -  s2(ti)] 

S2(t) 


,  t.  <  t  <_  t.+1,  i  =  0,  ....  m. 


In  this  formulation,  there  are  m  +  1  equally  spaced  points  in  time 
where  the  number  of  survivors  is  available. 

This  form  for  X2(t)  1S  convenient  when  life  table  data  are 
available.  This  may  often  be  the  case  since  many  experiments  have 
control  groups  of  the  species  under  investigation.  These  groups  are 
subjected  to  either  none  of  the  carcinogen  or  the  background  level  of  the 
carcinogen.  The  observed  hazard  rate  in  these  groups  will  be  a  good 
approximation  for  X^(t)  provided  that  the  background  contribution  of 
Aj(t;d)  is  small  as  will  often  be  the  case. 

The  specific  examples  used  throughout  this  dissertation  use  human 
life  table  data  to  obtain  X2(t).  The  data  are  from  Table  6.2  of  Vital 
Statistics  of  the  United  States  (National  Center  for  Health  Statistics 
(1984)).  The  contribution  of  X^(t;d)  to  these  mortality  data  is  assumed 
to  be  negl igible. 

The  cause-specific  probability  of  succumbing  to  cause  J  =  1  before 
time  t  is  given  by 

Pcs(t,l;d)  =  Jq  X^(u;d)  expt-A^u-.d)  -  A? ( u )  1  du  , 
or  equivalently,  in  a  more  convenient  computational  form. 


R? 


ifHiifci  '•J'JI'wl  kJ  x 


pcs(t*i;d)  = 


.£  {S2(ti)  expC-A1(t i ; d>TI  -  S2(ti+1)  exp[-A1(ti+1;d)] 


+  CS2(ti+i)  -  S2(ti)]  /t^+1  exp[-A^(u;d)]  du }  , 

with  t  =  tm+^.  The  computational  form  is  obtained  by  substituting  the 
piecewise  continuous  form  for  exp  [-A2(t)]  from  (4.1)  and  integrating  by 
parts  to  account  for  \^(t;d).  Depending  upon  the  parametric  form  for 
Aj(t;d),  the  integration  in  the  computational  form  may  be  done  directly, 
or,  a  numerical  integration  routine  may  be  used. 


»V*  fa’ .1  .*  **  , 


5.  RISK  CHARACTERIZATIONS  INVOLVING  TIME  TO  RESPONSE 


In  this  section  risk  characterizations  obtained  from  a  simple  model 
will  be  compared  to  improved  characterizations  obtained  from  the  cause- 
specific  model.  Two  risk  characterizations  involving  the  time  to 
response  are  considered;  namely,  the  probability  that  the  time  to  the 
response  precedes  a  specified  time  and  the  mean  response  free  period  in  a 
specified  period.  The  comparisons  will  be  made  in  terms  of  relative 
bi  as . 


5.1  Relative  Bias 

Relative  bias  refers  to  two  deterministic  representations  of  a 
non-random  risk  characteristic,  one  of  which  is  assumed  to  be  the  true 
representation . 

Definition  5.1  Let  and  ^  be  tw0  representations  of  a  risk 
characteristic  and  let  $2  136  true  representation.  The  relative 
bias  of  with  respect  to  $2  is 

RB(«1  ,®2)  =  ($1  “  =  V®2  "  1  • 

This  is  the  fraction  by  which  overstates  or  understates  Often 
lower  and  upper  bounds  can  be  calculated  for  the  relative  biases. 


5.2  Time-to-Response  Probability 


5.2.1  Bounds  on  Relative  Bias 


If  the  risk  characteristic  is  the  probability  that  the  time  to  a 
specified  response  precedes  a  specified  time  t  at  a  dose  d,  then  the 
relative  bias  of  P$(t;d)  with  respect  to  Pcs(t,l;d)  is 


RB(PS.PCS)  = 


p^TOld) 


In  this  case,  upper  and  lower  bounds  on  the  relative  bias  may  be 
calculated. 


Theorem  5.1  If  A^(t;d)  and  A2 ( t )  are  nondecreasing,  nonnegative, 
continuous  cumulative  hazard  functions  such  that  A^O-.d)  =  A?(0)  =  0, 
then  for  all  t  _>  0 

0  £  RB(PS,PCS)  £  exPCA2^)^  "  1  • 


Proof.  Since  A2(t)  is  nonnegative, 

exp[-A2(t)]  <_  1  for  all  t  >_  0  . 

Therefore, 

P s ( t ; d )  =  Jg  A1(u;d)exp[-A1(u;d)l  du 

-  /o  xi^u'»d)exPC-A1(,i;f1)]exp[-A2(u)]  du 
=  PrJt,l;d)  . 
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This  implies  RB(PS,PC$)  0  . 


Also,  the  hypothesis  on  A^t)  implies  that  expC-A^t)]  is  non¬ 
increasing  in  t.  Therefore, 

Pcs(t,l;d)  =  Jq  \1(u;d)exp[-A1(u;d)]exp[-A2(u)]  du 
>_  exp[-A2(t)]  /g  X^(u;d)exp[-A^(u ;d)]  du 


=  exp[-A2(t)]  P $ ( t ; d ) 


This  implies  that 


Ps(t;d) 


(tXdTlexpU2(t)] 


and  the  result  follows. 


The  conclusion  from  Theorem  5.1  is  that  the  simple  model  inflates 
the  probability  of  the  time  to  the  specified  response  preceding  t 
relative  to  the  cause-specific  probability.  This  is  true  not  only  for 
the  actual  probabilities  but  also  for  the  estimated  probabilities.  This 
is  because  A^(t;d)  from  the  simple  model  is  the  same  as  A^(t;d)  from  the 
cause-specific  model  as  shown  in  Section  3.  Thus,  the  proof  is  the  same 
with  the  true  cumulative  hazard  functions  replaced  by  the  estimated 
cumulative  hazard  functions. 


5.2.2  Attaining  the  Upper  Bound 

The  attainabi 1 ity  of  the  upper  bound  on  the  relative  bias  will  be 
discussed  in  terms  of  A^(t;d).  That  is,  A2(t)  is  considered  to  be  a 


19 


fixed  function  of  t  and  the  goal  is  to  find  Aj(t;d)  such  that  the  upper 
bound  in  Theorem  5.1  is  attained.  The  case  where  t  =  0  is  trivial  since 
the  lower  and  upper  bounds  are  equal  at  t  =  0.  Is  there  a  cumulative 
hazard  function  A^(t;d)  such  that  the  upper  bound  is  attainable  for  some 
t  >  0? 

Consider  the  case  where  the  carcinogen  has  no  effect  until  time  tg. 
Here,  tQ  is  called  the  latency  period.  In  this  case,  let  the  cumulative 
hazard  function  have  the  form 


Aj(t;d)  = 


P(t-tg) kg(d)  for  t  _>  t( 


otherwi se 


Then  let 


Fx(t  ;d)  = 


1  -  exp[-P (t-tg)Kg(d) ]  for  t  >_  tQ 


otherwise  . 


In  this  formulation,  as  P  increases,  the  carcinogen's  lethality 
increases.  For  all  t  tg, 

Pcs(t,l;d)  =  Jg  X1(u;d)exp[-A1(u;d)]exp[-A2(u)T  du 

=  Jg  exp[-A2(u)]  dF^u-.d)  . 


1$  dF^u-.d) 


RB(Ps,Pc$)  =  - - - i - 

/ g  exp[-A„(u)]  dF.(u;d) 


-  1  . 
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Thus,  as  g  ■*  ®, 


F  x  ( u ;  d ) 


1  for  u  >_  tg 


0  otherwise  , 


the  cdf  of  a  point  mass  at  tQ,  and 


/g  exp[-A2(u)]dF1(u;d)  >  exp[-A2(tQ)] 


This  implies 


RB(PS,PCS)  exp[A2(tQ)]  -  1  as  3 


Thus,  the  relative  bias  approaches  the  upper  bound  at  t  =  tn  as  M  °°  . 


This  example  shows  that  for  a  special  form  of  A.(t;d)  the  relative 


bias  approaches  the  upper  bound  at  some  t  >  0  in  the  limiting  case. 


However,  if  the  hypothesis  of  Theorem  5.1  is  changed  so  that  A2(t)  is 
additionally  required  to  be  strictly  increasing  then  a  strict  inequality 


holds  for  the  upper  bound.  That  is 


RB(Ps,PCs)  <  exp[A2(t)]  -  1 


5.2.3  Effect  of  Latency  Period 


Latency  period  affects  relative  bias  also.  To  investigate  this 


effect,  relative  biases  were  calculated  for  several  latency  periods.  The 


models  used  were 


P_  ( t ; d )  =  /*  \,(u;d)  exp[-A. (u ;d) ]  du 

5  1  1 
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Pcs(t,l;d)  =  J*  X.1(u;d)exp[-A1(u;d)  -  A2(u)]  du 
1/ 

where  A^t-.d)  =  p(t-tg)  g(d)  for  an  unspecified  g(d)  and  tQ  is  the 
latency  period.  For  each  t  =  60,  65,  70,  75,  80,  85,  90,  95,  100,  all 
combinations  of  the  parameter  k  and  latency  period  tg  were  considered 
where 

k  =  3,  4,  5,  6 
and 

tg  =  0,  10,  20,  30,  40,  50,  60  . 

These  ranges  for  the  parameters  k  and  tg  correspond  to  observed  estimates 
of  time-to-response  effects. 

The  life  table  data  were  truncated  at  t  =  85,  so  the  values  of  ^(t) 

are  extrapolated  from  t  =  85  to  t  =  100.  The  extrapolation  was  uone  by 

£ 

fitting  a  Weibull  cumulative  hazard  function,  A2(t)  =  yt  ,  to  the  life 
table  data.  A  least  squares  fit  gives  y  =  0.337Rxl0-11  so  that  for 
t  >  85,  A2(t)  =  yt6. 

The  values  of  Pxg(d)  were  chosen  such  that  Pcg(75,l;d)  =  10'6. 

This  helps  to  make  the  relative  biases  more  comparable  across  latency 
periods.  Therefore,  each  tg,  k  combination  uses  a  different  value  of 
8*g( d ) . 

Tables  3  through  6  give  the  relative  biases.  For  example  when 
k  =  3,  tg  =  60  and  t  =  100  the  relative  bias  of  the  simple  model 
representation  for  the  time-to-response  probability  with  respect  to  the 
cause-specific  model  representation  is  4.1023  which  implies  that 
P$(100;d)  is  5.1023  times  greater  than  P  ( 100 , 1 ; d ) ,  i.e.,  5.1023  times 
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Table  3 

Relative  bias  of  time-to-response  probability  for 
various  latency  periods  with  k  =  3 


upper 


t 

0 

10 

20 

30 

40 

50 

60 

bound 

60 

0.0901 

0.0987 

0.1093 

0.1227 

0.1405 

0.1650 

0.1990 

65 

0.1206 

0.1327 

0.1477 

0.1668 

0.1922 

0.2268 

0.2735 

0.3062 

70 

0.1647 

0.1820 

0.2035 

0.2310 

0.2671 

0.3159 

0.3829 

0.4742 

75 

0.2278 

0.2527 

0.2837 

0.3234 

0.3756 

0.4461 

0.5431 

0.7719 

80 

0.3198 

0.3566 

0.4028 

0.4623 

0.5414 

0.6497 

0.8029 

1.3875 

85 

0.4525 

0.5080 

0.5785 

0.6708 

0.7957 

0.9712 

1.2288 

2.7151 

90 

0.6275 

0.7098 

0.8157 

0.9565 

1.1505 

1.4296 

1.8504 

5.0207 

95 

0.8498 

0.9694 

1.1258 

1.3374 

1.6358 

2.0776 

2.7687 

10.9788 

100 

1.1210 

1.2907 

1.5161 

1.8276 

2.2787 

2.9698 

4.1023 

28.3121 

Table  4 

Relative 

bias  of 

time-to- 

response 

probability  for 

various 

latency  periods  with  k  = 

4 

t0 

upper 

t 

0 

10 

20 

30 

40 

50 

60 

bound 

60 

0.1017 

0.1103 

0.1206 

0.1334 

0.1497 

0.1710 

0.1990 

65 

0.1387 

0.1511 

0.1662 

0.1850 

0.2089 

0.2400 

0.2812 

0.3062 

70 

0.1931 

0.2113 

0.2335 

0.2611 

0.2960 

0.3410 

0.3992 

0.4742 

75 

0.2722 

0.2992 

0.3322 

0.3732 

0.4252 

0.4923 

0.5798 

0.7719 

80 

0.3910 

0.4324 

0.4835 

0.5477 

0.6301 

0.7385 

0.8849 

1.3875 

85 

0.5687 

0.6344 

0.7166 

0.8218 

0.9601 

1.1473 

1.4100 

2.7151 

90 

0.8123 

0.9149 

1.0451 

1.2150 

1.4431 

1.7592 

2.2119 

5.0207 

95 

1.1368 

1.2947 

1.4991 

1.7719 

2.1482 

2.6872 

3.4888 

10.9788 

100 

1.5542 

1.7926 

2.1082 

2.5407 

3.1577 

4.0797 

5.5266 

28.3121 
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Table  5 

Relative  bias  of  time-to-response  probability  for 
various  latency  periods  with  k  =  5 


20  30 


50  60 


0.1109 

0.1534 

0.2166 

0.3099 

0.4533 

0.6749 

0.9891 

1.4259 

2.0164 


0.1193 

0.1659 

0.2352 

0.3381 

0.4979 

0.7486 

1.1092 

1.6204 

2.3275 


0.1293 

0.1808 

0.2575 

0.3719 

0.5520 

0.8395 

1.2598 

1.8697 

2.7369 


0.1414 

0.1989 

0.2846 

0.4130 

0.6185 

0.9536 

1.4528 

2.1973 

3.2921 


0.1564 

0.2213 

0.3178 

0.4636 

0.7018 

1.1002 

1.7057 

2.6395 

4.0709 


0.1751 

0.2498 

0.3591 

0.5267 

0.8080 

1.2930 

2.0453 

3.2531 

5.2036 


0.2852 

0.4105 

0.6061 

0.9466 

1.5536 

2.5098 

4.1241 

6.9065 


upper 

bound 


0.1990 

0.3062 

0.4742 

0.7719 

1.3875 

2.7151 

5.0207 

10.9788 

28.3121 


Table  6 

Relative  bias  of  time-to-response  probability  for 
various  latency  periods  with  k  =  6 
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greater  than  the  actual  probability  of  observing  the  specified  response 
when  there  are  competing  risks.  Tables  3  through  6  show  the  increase  in 
relative  bias  as  the  latency  period  increases,  as  k  increases,  and  as  t 
increases.  For  example,  the  relative  bias  for  tg  =  40  and  k  =  5 
increases  from  0.1564  for  t  =  60  to  4.0709  for  t  =  100.  For  fixed  tg  and 
t,  say  40  and  75  respectively,  the  relative  bias  increases  from  0.3756 
for  k  =  3  to  0.4942  for  k  =  6.  Similarly,  for  k  =  5  and  t  =  75  the 
relative  bias  increases  from  0.3099  for  tg  =  0  to  0.6061  for  tg  =  60  . 

5.2.4  Fitting  a  Simple  Model  Without  a  Latency  Period  to  Data  from  a 
Cause-Specific  Model  with  a  Latency  Period 

The  relative  biases  given  previously  were  calculated  using  the 
simple  hazard  function  parameters.  Since  the  factor  L-^e.t)  in  the 
likelihood  associated  with  the  cause-speci f_ic  model  has  the  same  form  as 
the  likelihood  associated  with  the  simple  model,  the  maximum  likelihood 
estimates,  5,  are  the  same  under  both  models.  Hence  Tables  3  through  6 
also  give  the  relative  bias  of  the  estimated  simple  model  relative  to  the 
estimated  cause-specific  model.  However,  if  data  are  from  a 
cause-specific  model  with  the  cumulative  hazard  function  for  the 
specified  response  given  by 

P(t-tg)  kg(d)  for  t  >_  tg 

A1(t;d)  =< 

0  otherwise  , 

and  these  data  are  fit  by  a  simple  model  with 


»] 

ft 


A1(t;d)  =  0  t*  g(d) 


then  the  likelihoods  do  not  have  the  same  form. 

A  small  Monte  Carlo  study  was  conducted  to  obtain  estimates  of  the 
relative  biases  associated  with  not  only  using  simple  models  but  also 
omitting  the  latency  period.  In  this  study  the  probability  of  the 
specified  response  occurring  before  time  t  would  be  of  the  order  10”^ 
if  the  dose  levels  of  primary  concern  were  used.  This  would  require  that 
the  number  of  simulations  be  very  large  in  order  to  insure  that  at  least 
a  few  specified  responses  occurred.  To  alleviate  this  problem,  the  study 
was  conducted  using  doses  high  enough  to  give  moderate  to  large 
percentages  of  the  outcomes  being  the  specified  response.  The  estimated 
probabilities  were  then  extrapolated  down  to  the  low  dose  region.  The 
simulation  procedure  is  as  follows: 

The  experimental  data  were  generated  using 

?cs(t,l:d)  =  /£  ^(ujdJexpC-Ajfujd)  -  A2(u)]  du 


Aj(t;d)  =  8(t-tgj  g(d)  ,  t  >  tQ  and  tQ  =  40 


Also,  let 


Ps(t;d)  =  fg  \  (u;d)exp[-Aj(u;d)]  du 


A*  (t;d)  =  6*tk  g(d)  . 


with 
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The  dose  function,  g(d),  was  not  specified;  however,  it  was  assumed 
that  at  the  dose  level  d^  where  the  estimated  relative  bias  was  to  be 
computed,  g ( dQ )  =  1  and  pcs(75,l;dg)  =  10~6.  This  implied  a  value  for 
p.  The  experimental  data  were  simulated  with  1000  observations  at  each 


dj,  ...,  dg  where 


Pcs^max**:t,i^  “  ^i 


(5.1) 


with  fj  =  0.3,  f2  =  0.4,  fj  =  0.5,  f4  =  0.6,  fg  =  0.7  and  tmax  =  125 
corresponding  to  a  maximum  possible  lifetime.  Equation  (5.1)  implied  the 
val ue  of  g(d^ ) . 

For  each  of  300  sets  of  experimental  data,  the  maximum  likelihood 

A*  A* 

estimates  p  ,  and  k  were  determined  as  well  as  RR(ps»pcs)*  Tables  7 
through  10  give  the  sample  means  of  the  RR(^S »pcs) 's.  The  variances  are 
not  given  in  the  tables  but  they  were  small.  Several  95%  confidence 
intervals  were  constructed  and  resulted  in  interval  half-widths  which 
were  about  2%  of  the  relative  bias  estimates.  The  last  column  in  each 
of  these  tables  is  the  corresponding  column  from  Tables  3  through  6  with 
tQ  =  40.  This  last  column  shows  the  actual  relative  bias  of  the  simple 
model  with  latency  period  relative  to  the  cause-specific  model. 

Tables  7  through  10  show  that  the  effect  of  not  including  a  latency 
period  in  a  simple  model  when  the  underlying  cause-specific  model  has  a 
latency  period  usually  increases  the  relative  bias,  especially  for  larger 
k.  The  difference  is  not  so  great  for  the  middle  values  of  t,  say  t  =  75 
or  80,  as  it  is  for  smaller  or  larger  values  of  t.  For  example,  for 
k  =  5  and  fg  =  0.5  the  relative  bias  without  a  latency  period  is  0.3684 


©W 3SS 


j- 
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Table  7 

Relative  bias  of  time-to- response  probability 
for  k  =  3  and  tg  =  40. 


Known  relative 
bias  of  the 

Estimated  relative  bias  of  the  simple  model  simple  model 

without  a  latency  period  with  a 

latency  period 


t 

fl  =  0.3 

f2  =  0.4 

f3  =  0.5 

VO 

• 

o 

II 

4- 

f5  =  0.7 

60 

0.2755 

0.2257 

0.1782 

0.1197 

0.0671 

0.1405 

65 

0.1553 

0.1263 

0.1035 

0.7447 

0.0616 

0.1922 

70 

0.1570 

0.1429 

0.1397 

0.1350 

0.1595 

0.2671 

75 

0.2458 

0.2457 

0.2628 

0.2842 

0.3533 

0.3756 

80 

0.4306 

0.4467 

0.4894 

0.5446 

0.6756 

0.5414 

85 

0.7456 

0.7840 

0.8634 

0.9683 

1.1941 

0.7957 

90 

1.2218 

1.2932 

1.4282 

1.6098 

1.9848 

1.1505 

95 

1.9231 

2.0453 

2.2663 

2.5687 

3.1816 

1.6358 

100 

2.9295 

3.1299 

3.4842 

3.9761 

4.9663 

2.2787 

Table 

8 

Relative  bias  of  time-to- response 
for  k  =  4  and  tg  =  40. 

probability 

t 

Estimated  relative 
without 

fx  =  0.3  f2  =  0.4 

i  bias  of  the  simple 
a  latency  period 

fg  =  0.5  f^  =  0.6 

model 

f5  =  0.7 

Known  relative 
bias  of  the 
simple  model 
with  a 

latency  period 

60 

0.5602 

0.4642 

0.3714 

0.2663 

0.1666 

0.1497 

65 

0.3068 

0.2515 

0.2034 

0.1491 

0.1132 

0.2089 

70 

0.2515 

0.2211 

0.2031 

0.1850 

0.2025 

0.2960 

75 

0.3194 

0.3097 

0.3200 

0.3382 

0.4179 

0.4252 

80 

0.5154 

0.5276 

0.5727 

0.6381 

0.8071 

0.6301 

85 

0.8848 

0.9297 

1.0267 

1.1652 

1.4807 

0.9601 

90 

1.4835 

1.5784 

1.7595 

2.0193 

2.5848 

1.4431 

95 

2.4334 

2.6116 

2.9347 

3.4034 

4.4075 

2.1482 

100 

3.9150 

4.2345 

4.8000 

5.6312 

7.4077 

3.1577 
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Table  9 

Relative  bias  of  time-to-response  probability 
for  k  =  5  and  tg  =  40. 

Known  relative 
bias  of  the 


Estimated  relative  bias  of  the  simple  model  simple  model 

without  a  latency  period  with  a 

latency  period 


t 

f x  =  0.3 

f2  -  0.4 

f3  =  0.5 

VO 

• 

o 

II 

4- 

f5  =  0.7 

60 

0.8866 

0.7307 

0.5867 

0.4250 

0.2572 

0.1564 

65 

0.4644 

0.3818 

0.3047 

0.2236 

0.1556 

0.2213 

70 

0.3428 

0.3002 

0.2616 

0.2316 

0.2388 

0.3178 

75 

0.3867 

0.3750 

0.3684 

0.3866 

0.4792 

0.4636 

80 

0.5899 

0.6117 

0.6423 

0.7230 

0.9422 

0.7018 

85 

1.0133 

1.0832 

1.1704 

1.3525 

1.7929 

1.1002 

90 

1.7357 

1.8855 

2.0696 

2.4312 

3.2775 

1.7057 

95 

2.9572 

3.2497 

3.6106 

4.3057 

5.9278 

2.6395 

100 

5.0116 

5.5661 

6.2578 

7.5852 

10.7137 

4.0709 

Table 

10 

Relative  bias  of  time-to-response 
for  k  =  6  and  tg  =  40, 

probability 

t 

Estimated  relative 
without 

fl  =  0.3  f2  =  0.4 

i  bias  of  the  simple 
a  latency  period 

f3  =  0.5  f4  «  0.6 

model 

f5  -  0.7 

Known  relative 
bias  of  the 
simple  model 
with  a 

latency  period 

60 

1.2522 

1.0328 

0.7871 

0.5755 

0.3529 

0.1614 

65 

0.6233 

0.5105 

0.3856 

0.2851 

0.1950 

0.2308 

70 

0.4269 

0.3656 

0.3029 

0.2665 

0.2677 

0.3349 

75 

0.4436 

0.4180 

0.4032 

0.4250 

0.5275 

0.4942 

80 

0.6513 

0.6619 

0.7017 

0.8002 

1.0575 

0.7609 

85 

1.1220 

1.1847 

1.3100 

1.5391 

2.0820 

1.2211 

90 

1.9608 

2.1141 

2.3937 

2.8675 

3.9685 

1.9411 

95 

3.4585 

3.7847 

4.3654 

5.3269 

7.5760 

3.1039 

100 

6.1521 

6.8226 

8.0128 

9.9777 

14.6566 

4.9969 
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for  t  =  75  (compared  to  0.4636  with  a  latency  period),  but  0.5878  for 
t  =  60  compared  to  0.1564)  and  6.2578  for  t  =  100  (compared  to  4.0709). 
The  relative  biases  for  the  simple  model  without  a  latency  period  tend  to 
decrease  as  the  values  of  f-  increase  for  smaller  t  but  the  relative 
biases  increase  as  the  values  of  f ^  increase  for  larger  t. 

5.3  Mean  Response  Free  Period 
5.3.1  New  Definitions 

Mean  response  free  period  was  introduced  as  a  risk  characterization 
by  the  Society  of  Toxicology  EDg^  Task  Force  (1981).  It  is  defined  as 
follows: 

Definition  5.2  The  mean  response  free  period,  MFP(t;d),  is  the 
expected  time  at  dose  d  without  a  specified  response  in  the  first  t 
units  of  time.  Analytically, 

MFP(t;d)  »  t  -  j|j  (t-u)  [£i(Jjil]T=u  du  . 

Using  the  simple  model,  this  becomes 

MFPs(t;d)  =  t  -  /g  (t-u)\^(u;d)exp[-A^(u;d)]  du  .  (5.2) 

The  concept  of  a  mean  response  free  period  must  be  refined  to 
incorporate  the  multiple  possible  responses  in  the  cause-specific  model. 
Definition  5.3  The  total  mean  response  free  period  under  the 
cause-specific  model  is 


MFPcs(t;d)  =  t  -  /q  (t-u)(X1(u;d)+X2(u) )exp[-A1(u;d)-A2(u)]  du. 

This  is  the  expected  time  without  a  response  due  to  cause  1  or  2  in 
the  first  t  units  of  time  at  dose  d.  It  is  also  possible  to  define  mean 
response  free  periods  associated  with  each  cause. 

Definition  5.4  The  cause-specific  mean  response  free  period 
associated  with  cause  j  is  the  expected  time  without  a  response  due 
to  cause  j  in  the  first  t  units  of  time  at  dose  d.  Analytically, 

MFPcs (t ,j ;d)  =  t  -  /q  (t-u)Xj.(u;d)exp[-A1(u;d)-A2(u)]  du  . 

The  relationship  among  the  mean  response  free  periods  under  the 
cause-specific  model  is 

MFPftjd)  =  MFP  (t,l;d)  +  MFP  ft,2;d)  -  t  . 


Since  the  simple  model  attempts  to  model  the  probability  of  the 
specified  response,  comparisons  will  be  made  between  MFP$(t;d)  and 
MFP  (t,l;d)  and  measured  in  terms  of  relative  bias. 


5.3.2  Bounds  on  Relative  Bias 


From  Definition  5.1,  the  relative  bias  of  the  simple  mean  response 
free  period,  equation  (5.2),  with  respect  to  the  cause-specific  mean 
response  free  period  due  to  cause  1,  Definition  5.4,  is 

,  ,  MFP  ( t ; d ) 

R8(MFP5(t;d,,MFPc5(t.l;d))  =  _L__  -  1  . 


In  Section  5.2.1,  it  was  shown  that  the  simple  model  overstates  the 
probability  of  a  specified  response  before  time  t.  Therefore,  it  follows 
that  the  simple  model  should  understate  the  mean  response  free  period 
associated  with  cause  1.  The  following  theorem  verifies  this  conjecture. 

Theorem  5.2  If  A^tjd)  and  A2(t)  are  nondecreasing,  nonnegative 
cumulative  hazard  functions  such  that  A^(0;d)  =  A^fO)  =  0.  Then  for  all 


t  >  0, 


-1  <  RB(MFP  (t;d) ,MFP  (t,l;d) )  <  0 


Proof:  Since  A2U)  is  nondecreasing  in  t, 
/q  (t-u)A1(u;d)exp[-A1(u;d)]  du 


i  Jg  (t-u)X1(u;d)exp[-A1(u;d)lexp[-A2(u)]  du  , 


so  that 


MFPs(t;d)  <  MFPc$(t,l;d)  . 


This  implies 


RB(MFPs(t;d),MFPcs(t,l;d))  <  0 


Also,  MFPcs(t,l;d)  may  be  written  as 


MFPcs(t,l;d)  =  t[l  -  Pcs(t,l;d)] 


+  /g  ux^j(u;d)exp[-A1(u;d)  -  A2(u)]  du 
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since  both  terms  are  nonnegative.  Similarly  MFP$(t;d)  0.  Therefore, 
RB(MFPs(t;d),MFPcs(t,l;d))  >  -1  . 

5.3.3  Examples 

Some  examples  are  presented  to  illustrate  the  magnitude  of  the 
relative  bias  of  the  simple  mean  response  free  period  with  respect  to  the 
total  mean  response  free  period  under  the  cause-specific  model.  Using 
the  hazard  function  for  the  competing  risk,  described  in  Section  4, 
computational  forms  for  MFP  { t ;d) ,  MFP  (t;d),  MFP  (t,l;d),  MFP  (t ,2;d) 
and  MFP$(t;d)  can  be  obtained  in  a  straightforward  manner  sometimes 
including  one  or  two  integrations-by-parts.  They  are 

m 

MFPcs(t;d)  «  t  -  t  l  (S2(ti )expC-A1(ti Jd)3 

-  S2(t.+1)exp[-A1(ti+1;d)]} 
m 

+  l  (tiS2(ti)exp[-A.(ti ;d)] 
i=0  1  c  1  11 

-  ti+1S2(ti+1)exp[-A1(ti+1;d)]} 

m  t.  . 

+  l  (a.  L  1  u*exp[-A. (u ;d) 1  du} 
i=0  1  li  1 

m  t-., 

+  l  (b.  f.1  1  exp[-A.(u;d)]  du)  , 
i=0  li  1 


MFP  (t,2;d)  =  t  +  t  l  (a.  /.  exp[-A. ( u ;d ) ]  du) 

cs  1-0  1  ci  1 


-  I  (a,  L1  +  1  uxexp[-A. (u;d)]  du}  , 
i-0  i 


MFPcs(t,l:d)  =  MFPc$(t;d)  -  MFPcs(t,2;d)  +  t  , 


MFP$(t;d)  =  txexp[-A^(t;d)l 


+  j  expC-Aj (t. ;d)1  -  t^+1expC-Aj(t^+j ;d)l 


/ti+1  exp[-A^(u;d)]  du} 


with  t  =  t 


These  formulae  are  used  to  obtain  the  mean  response  free  periods 
used  to  calculate  RB(MFPs(t ;d) ,MFPcs(t , 1 ;d) )  in  Tables  11  through  14. 

The  form  used  for  the  cumulative  hazard  function  due  to  the  cause  of 
interest  was  the  same  as  in  Section  5.2,  i.e.,  A^(t;d)  =  pt  g(d). 

Tables  11  through  14  correspond  to  k  =  3,  4,  5  and  6,  respectively.  Eac 
column  in  these  tables  is  calculated  for  a  different  Pxg(d)  where  Pxg(d) 
is  chosen  so  that 


Pcs(75,l;d)  =  fQ  (  =  lO-0) 


Pcs^tmax,1;(1^  f  i  ’  1  =  1 . 


where  f^  =  0.3.  f„  =  0.4.  f =  0.5.  f„  =  0.6.  fc  =  0.7  and  t _  =  125, 


Table  12 
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The  resulting  mean  response  free  periods  are  not  given  in  the  tables 
but  both  the  simple  and  cause-specific  mean  response  free  periods 
associated  with  fg  were  very  close  to  the  value  for  t.  However,  the  mean 
response  free  periods  decreased  as  f.  increased.  For  example, 
MFPcs(100;d)  =  99.9999  and  MFPs(100;d)  =  99.9999  for  fQ  =  lCf6  and  k  =  3; 
but  MFPcs(100;d)  =  78.3223  and  MFPs(100;d)  =  72.6925  for  f3  =  0.5  and 
k  =  3. 

Some  of  the  conclusions  that  can  be  drawn  from  the  tables  are  as 
follows: 

(i)  The  magnitude  (absolute  value)  of  the  relative  bias  increases 
as  t  increases. 

(ii)  For  smaller  values  of  t,  the  magnitude  of  the  relative  bias 
decreases  as  k  increases,  but  for  large  t  the  magnitude 
increases  as  k  increases. 

(iii)  The  magnitude  of  the  relative  bias  tends  increases  as  the 

probability  of  a  specified  response  over  a  lifetime  increases. 

In  comparison  to  the  relative  biases  in  Tables  3  through  6  for  the 
time-to-response  probabilty,  the  relative  biases  in  Tables  11  through  14 
for  the  mean  response  free  period  are  small.  For  example,  if  the 
probability  of  a  specified  response  before  time  t  =  75  is  10  ^  then  the 
relative  bias  in  Table  4  for  the  time-to-response  probability  for  k  =  4, 
tg  =  40  and  t  =  75  is  0.4252;  but  the  relative  bias  in  Table  12  for  the 
mean  response  free  period  for  k  =  4,  and  t  =  75  is  only  -0.3146xl0~'7. 
Thus,  P s { t ; d)  is  1.4252  times  greater  than  P cs ( t ; d )  but  the  mean  response 
free  periods  are  almost  the  same  under  both  models. 
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6.  RISK  CHARACTERIZATIONS  INVOLVING  DOSE 


In  this  section,  as  in  Section  5,  risk  characterizations  obtained 
from  a  simple  model  will  be  compared  to  corresponding  risk 
characterizations  obtained  from  a  cause-specific  model.  Two  risk 
characterizations  involving  dose  are  considered;  namely,  mean  free  dose 
and  virtually  safe  dose.  The  comparisons  will  be  made  in  terms  of 
relative  bias. 


6.1  Mean  Free  Dose 


6.1.1  Definition 


The  mean  free  dose  was  introduced  by  the  Society  of  Toxicology  EDQ1 
Task  Force  (1981).  It  is  interpreted  as  the  dose  which  would  result  in  a 
specified  fractional  reduction  in  the  mean  response  free  period. 

Definition  6.1  The  mean  free  dose  at  time  t  is  the  dose  d  which 
satisfies 


MFP(t;d)  =  1  _ 
MFP(t ;0) 


with  0  <  e  <  1  and  e  being  near  0.  The  mean  free  dose  at  time  t  is 
denoted  by  MFD(t;e). 

This  definition  may  be  used  with  any  of  the  mean  response  free 
periods  defined  in  Section  5.3.  The  notation  for  a  mean  free  dose  will 
correspond  to  the  notation  used  in  Section  5.3  for  mean  response  free 


period.  For  example,  if  Definition  6.1  uses  MFP  (t,l;d),  then  the 


notation  for  the  corresponding  mean  free  dose  would  be  MFDcs(t,l;e) 


ilU 


6.1.2  Approximations  of  Mean  Free  Dose 

By  assuming  a  fairly  general  form  for  A^tjd),  simple  approximations 
are  available  for  MFD  (t,l;e)  and  MFD  (t;e).  They  can  be  derived  as 

v  j  b 

follows:  Let 

Aj(t;d)  =  (1  +  dr)  ptk  . 

Then  MFDcs(t,l;e)  is  the  dose  d  that  satisfies 

MFPcs(t,l;d)  =  MFPcs(t,l;0)  -  e*MFPcs(t ,1 ;0)  , 

with 

MFPcs(t,l;d)  = 

t  -  /£  (t-u)(l+dr)Bkuk“1exp[-(l+dr)puk  -  A2(u)]  du  .  (6.1) 

Using  the  approximation  for  small  p,  exp(-drpuk)  =  1,  equation  (6.1) 
becomes 

MFPc$(t,l;d)  = 

t  -  /*  (t-u) (l+dr)pkuk-1exp[-puk  -  A2(u)]  du  . 

1 

♦ 

Thus  MFDc$(t,l;e)  is  approximately  the  dose  d  which  satisfies  j 

i 

t  -  /J  (t-u) (l+dr)pkuk_1exp[-Puk  -  A2(u)]  du  j 

) 

i 

=  t  -  /q  (t-u)Bkuk_1exp[-puk  -  A2(u)]  du  ! 

I 

-  e*MFPcs(t,l;0)  ,  j 

I 

I 

I 

I 

V 
1 
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which  can  be  rewritten  as 


(l+dr)  / q  (t-u^ku^expC-pu1*  -  A£(u)]  du 
=  /«  (t-tOPku^expC-pu1*  -  a9(u)]  du 


+  e*MFPcs(t,l;0)  . 


Solving  for  the  dose  d  the  result  is  that 


.  e*MFP  (t.l;0) 

MFD  (t  ,1  ;e )  =  - - - 

CS  t-MFP  (t,l;0) 


(6.2) 


A  similar  derivation  for  MFD$(t;e)  yields 


exMFP  (tfO) 

MFD  ft;e)  =  - 1- - 

s  t-MFP  (t;0) 


(6.3) 


Tables  15  through  22  give  examples  of  the  mean  free  doses  calculated 
from  numerical  search  routines  and  the  approximations  (6.2)  and  (6.3). 
These  mean  free  doses  were  calculated  for  t  =  75  and  e's  that  correspond 
to  losses  of  1  hour,  1  day,  1  week  and  1  month  in  a  75  year  period.  The 
values  of  r  were  r  =  1  and  r  =  3  corresonding  to  linear  and  nonlinear 
dose  functions,  respectively. 

Generally  the  approximations  agree  with  the  mean  response  free  dose 


obtained  by  the  numerical  search  technique  to  about  three  digits. 
However  as  the  dose  increases  the  approximation  is  not  as  good.  For 
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Table  15 

Approximations  of  mean  free  dose  with  k  =  3  and  linear  dose  function 


MFDs(t;e)  MFDcs(t,l;e) 


Using 

Using 

Search 

Search 

e 

Technique 

Approximation 

Technique 

Approximati on 

1.52*10*5 

0.4955*10+1 

0.4955*10+1 

0.5507 *10+1 

0.5505*10+1 

3.65*10*^ 

0.1189*10+3 

0.1189*10+3 

0.1322*10+3 

0.1322*10+3 

2.55*10"4 

0.8327*10+3 

0.8325*10+-3 

0.9255*10+3 

0.9252*10+3 

1.11*10'2 

0.3624*10+4 

0. 3620*10+4 

0.4028*10+4 

0.4023*10+4 

Table  16 

Approximations  of  mean  free  dose  with  k  =  3  and  nonlinear  dose  function 


1 


1 

r 

MFD 

$  (t  »e ) 

MFDcs(t,l;e)  1 

► 

r 

Using 

Using 

t 

Search 

Search 

£ 

Technique 

Approximati on 

Technique  Approximation 
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Table  17 

Approximations  of  mean  free  dose  with  k  =  4  and  linear  dose  function 


MFDs(t;e) 


MFDcs(t,l;e) 


'4 

Using 

Using 

Search 

Search 

4 

. 

e 

Technique 

Approximation 

Technique 

Approximation 

1.52x10"6 

0.5978x10+1 

0.5978x10+1 

0.6821x10+1 

0.6821 x10+1 

3.65*10“5 

0. 1435x10+3 

0. 1435x10+3 

0.1637x10+3 

0. 1637x10+3 

2.55x10“4 

0. 1005x10+4 

0.1004x10+4 

0 . 1146  x10+4 

0.1146x10+4 

l.llx10‘2 

0.4374x10+4 

0.4367x10+4 

0.4991 x10+4 

0.4983x10+4 

Table  18 

Approximations  of  mean  free  dose  with  k  =  4  and  nonlinear  dose  function 


MFD  ft;e) 


MFD(t,l;e) 


1 

4 

Us  i  ng 

Using 

4 

Search 

Search 

k  e 

V 

Technique 

Approximati on 

Techni que 

Approximation 

Table  19 


Approximations  of  mean  free  dose  with  k  =  5  and  linear  dose  function 


MFD  (t;e) 


Search 

Technique  Approximation 


MFD_(t,l;e) 


Usi 

Search 
Technique 


Approximati on 


1.52x10"6 

0.6967X10+1 

0.6967xlO+1 

0.8144x10+1 

0.8144x10+1 

3.65*10-5 

0.1672x10+3 

0. 1672x10+3 

0.1955x10+3 

0.1955x10+3 

2.55x10‘4 

0 . 1 171 x 10+4 

0.1171x10+4 

0.1369x10+4 

0.1368x10+4 

l.llxlO'2 

0.5099*10+4 

0.5090x10+4 

0.5961x10+4 

0.5949x10+4 

Table  20 

Approximations  of  mean  free  dose  with  k  =  5  and  nonlinear  dose  function 


MFD_ (t ;e ) 


Search 

Technique  Approximation 


MFD  (t ,  1 ;  e  ) 


Using 

Search 

Technique  Approximation 


1.52x10 

0.1910x10  1 

0.1910x10 

0.2012x10  1 

0.2012x10 

3.65x10-5 

0.5509x10+1 

0.5509x10+1 

0.5804x10+1 

0.5803x10 

2.55x10-4 

0. 1054x10+2 

0.1054xl0+2 

O.lllOxlO*2 

O.lllOxio 

l.llx10'2 

0. 1721x10+2 

0.1720xl0+2 

0. 1831 x10+2 

0.1812x10 

'V 
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Table  21 

Approximations  of  mean  free  dose  with  k  =  6  and  linear  dose  function 


MFDs(t;e) 


Using 

Search 

Technique  Approximation 


MFDcs(t,l;e) 


Using 

Search 

Technique  Approximation 


1.52x10'° 

0.7932x10  1 

0.7932x10  A 

0.9478x10 

0.9478x10  1 

3.65xl0'5 

0.1904xl0+3 

0.1904xl0+3 

0.2275xl0+3 

0.2275xl0+3 

2.55xl0'4 

0.1333xl0+4 

0.1333xl0+4 

0.1593xl0+4 

0.1592xl0+4 

l.llxlO'2 

0.5807xl0+4 

0.5795xl0+4 

0.6940xl0+4 

0.6924xl0+4 

Table  22 

Approximations  of  mean  free  dose  with  k  *  6  and  nonlinear  dose  function 


MFD(t;e) 


Using 

Search 

Technique  Approximation 


MFD„.(t,l;e) 


Using 

Search 

Technique  Approximation 


1.52x10'° 

0.1994x10 

0.1994x10 

0.2116x10 

0.2116x10 

3.65xl0"5 

0.5753xl0+1 

0.5753xl0+1 

0.6105xl0+1 

0.6104x10 

2.55xl0’4 

0.1101xl0+2 

0.1100xl0+2 

0.1168xl0+2 

0.1168x10 

l.llxlO'2 

0.1797xl0+2 

0.1796xl0+2 

0.1907  xl0+2 

0.1906x10 
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example  at  k  =  4  and  t  =  1.52*10"^  with  a  linear  dose  function  both  the 

approximation  and  the  search  method  yield  a  simple  mean  response  free 

dose  of  0.5978*10+*  but  at  e  =  1.11*10”^  the  approximation  is 
+4 

0. 4367*10  while  the  search  technique  gives  a  simple  mean  response  free 

dose  of  0.4374*10+^.  This  is  as  expected  since  the  approximation  of 
p  k  * 

exp(-d  Pu  )  =  1  is  better  for  low  doses. 

For  the  same  reason  the  approximation  is  not  as  good  as  k 
increases.  This  can  be  seen  by  comparing  the  approximation  and  search 

_p 

values  at  e  =  l.llxio  in  Table  15  with  the  corresponding  values  in 
Table  21. 

Also,  the  approximation  is  better  for  the  nonlinear  dose  function. 
This  is  because  the  nonlinear  dose  function  results  in  a  lower  mean 
response  free  dose  which  in  turn  makes  the  approximation  exp(-dr|3uk)  =  1 
better. 

6.1.3  Relative  Bias 

The  relative  bias  will  depend  upon  which  form  of  the  mean  free  dose 
is  considered.  Since  any  change  in  the  total  mean  response  free  period 
due  to  the  introduction  of  a  carcinogen  would  be  of  concern,  the  total 
mean  free  dose  is  explicitly  considered  in  this  subsection.  The  mean 
free  doses  are  calculated  from  the  simple  mean  response  free  period, 
equation  (5.2),  and  the  total  mean  response  free  period  under  the 
cause-specific  model.  Definition  5.3.  The  relative  bias  of  simple  mean 
free  dose  with  respect  to  total  mean  free  dose  under  the  cause-specific 


model  is 


RB(MFD,MFD  J  = 
s  cs 


MFDs(t;e) 

MFDcs(t;£) 


1  . 


The  form  of  the  specified  response's  cumulative  hazard  function  used 
to  calculate  the  mean  free  dose  is 


Aj(t;d)  =  (1  +  dr)(3tk  , 


(6.4) 


with  r  =  1  and  r  =  3  corresponding  to  linear  and  nonlinear  dose 
functions,  respectively.  The  value  of  p  is  such  that  P  ( 75, 1 ;0)  =  10“6, 
and  the  mean  free  doses  are  calculated  at  t  *  75.  Tables  23  and  24  give 
the  relative  biases  for  the  linear  and  nonlinear  dose  functions  with 
k  =  3,  4,  5  and  6.  The  same  values  of  e  as  those  in  Tables  15  through  22 
were  used  to  compute  the  relative  biases. 

Tables  23  and  24  show  the  following: 

(i)  All  relative  biases  are  negative,  indicating  that  MFDs (t ; e ) 
understates  MFDcs(t;s). 

(ii)  Relative  biases  are  virtually  unaffected  by  changes  in  e.  This 
is  expected  since  e  cancels  when  the  relative  bias  is 
calculated  using  the  approximation. 

(iii)  Relative  biases  have  greater  absolute  values  as  k  increases, 

(iv)  The  magnitude  of  the  relative  bias  is  greater  when  the  dose 
function  is  linear  than  when  it  is  nonlinear. 

The  effect  of  using  the  simple  model  instead  of  the  true  cause- 
specific  model  is  larger  when  the  dose  function  is  linear.  For  k  =  6  and 


a  linear  dose  function  the  relative  bias  is  -0.1925.  This  means  that  the 


Table  23 


Relative  bias  for  mean  free  dose  with  linear  dose  function 


k 


c 

3 

4 

5 

6 

1 .52x10”6 

-0.1301 

-0.1551 

-0.1755 

-0.1925 

3.65*10'5 

-0.1301 

-0.1551 

-0.1755 

-0.1925 

2.55*10‘4 

-0.1302 

-0.1552 

-0.1756 

-0.1925 

1.11*10"2 

-0.1302 

-0.1553 

-0.1757 

-0.1927 

Relative  bias 

e 

Table  24 

for  mean  free  dose  with  nonlinear  dose 

k 

3  4  5 

function 

6 

1.52*10-6 

-0.04541 

-0.05464 

-0.06231 

-0.06878 

3.65xl0'5 

-0.04541 

-0.05465 

-0.06231 

-0.06878 

2.55x10-4 

-0.04542 

-0.05465 

-0.06232 

-0.06880 

l.llxlO"2 

-0.04544 

-0.05469 

-0.06238 

-0.06887 

MFDs(t;e)  is  about  0.8  times  as  large  as  MFDcs(t;e).  That  is,  the  mean 
response  free  dose  calculated  under  the  simple  model  is  only  0.8  times  as 
large  as  the  true  mean  response  free  dose  when  competing  risks  are 
considered. 


v  -■ 


6.2  Virtually  Safe  Dose 


6.2.1  Definition 

The  virtually  safe  dose  was  considered  in  Hartley  and  Sielken  (1977) 
but  was  used  with  quantal  response  models  prior  to  1977.  It  is 
interpreted  as  the  dose  which  would  result  in  a  specified  slight  increase 
in  the  probability  of  the  specified  response  occurring  before  time  t  over 
this  probability  at  dose  0. 

Definition  6.2  The  virtually  safe  dose  at  time  t  is  the  dose  d 
that  satisfies 

P(t;d)  =  P(t;0)  +  n  , 

with  0  <  n  <  1  and  it  near  0.  The  virtually  safe  dose  is 
denoted  by  VSD(t;it). 

Since  the  probabilities  that  were  compared  under  the  simple  and 
cause-specific  models  were  P$(t;d)  and  Pcs(t,l;d),  in  Section  5,  the 
virtually  safe  doses  are  calculated  in  terms  of  these  two  probabilities 
using  Definition  6.2  and  the  corresponding  virtually  safe  doses  will  be 
VSD$(t;it)  and  VSDcs(t,l;it)  respecti vely.  The  VSDcs(t,l;u)  refers  to 
increasing  the  cause-specific  probability  of  the  specified  response  as 
opposed  to  the  total  probability  of  some  response. 

6.2.2  Approximations  of  Virtually  Safe  Dose 

As  with  the  mean  response  free  dose,  simple  approximations  for  the 
virtually  safe  doses  are  available  if  a  fairly  general  form  for  A.^(t;d) 


is  assumed.  Let 


I 


A1(t;d)  =  (1  +  dr)  f?th 


Then  VSDC$ ( t ,  1 ; it)  is  the  dose  d  that  satisfies 


/q  (l+dr)pkuk“1exp[-{l+dr)puk  -  a2(u)]  du 


=  / g  Pkuk~'''exp[-(3uk  -  A^u)]  du  +  it 


Using  the  approximation,  for  small  p,  exp(-drpu*)  =  1  this  becomes 


(l+dr)  /g  pkuk_1exp[-puk  -  A^(u )]  du 


=  Jg  Pkuk  ^exp[-puK  -  A2(u)]  du  +  it 


Solving  for  the  dose  d  the  result  is  that 


VSD  (t  ,1  ;it)  «  - - - 

cs  Pcs(t.i;0) 


A  similar  derivation  for  VSD$  yields 


VSD  (t  ;u)  =  - - - 

s  Ps(t;0) 


Tables  25  through  32  give  examples  of  the  virtually  safe  doses 

calculated  from  a  numerical  search  routine  and  the  approximations  (6.5) 

and  (6.6).  The  virtually  safe  doses  are  calculated  for  t  =  75,  n  =  10“ 
-6  -5  -4 

10  ,  10  and  10  .  The  values  of  r  were  the  same  as  for  the  mean 


response  free  dose;  namely,  r  =  1  and  r  =  3. 


V*  -> 


Table  25 

Approximations  of  virtually  safe  dose  with  k  =  3 
and  linear  dose  function 


VSDs(t;n) 


Using 

Search 

Technique  Approximation 


VSDcs(t,l;it) 


Us  i 
Search 

Technique  Approximation 


j 

Using 

Search 

'Jsi  ng 
Search 

1 

Technique 

Approximation 

Technique 

Approximati on 

I 
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Table  27 

Approximations  of  virtually  safe  dose  with  k  =  4 
and  linear  dose  function 


VSDs(t;*) 


VSDcs(t,l;«) 


u 

Using 

Search 

Technique 

Approximation 

Using 

Search 

Technique 

Approximation 

10"7 

0.786OX10'1 

0.7860X10"1 

0.1000*10+0 

0.1000xl0+0 

10-6 

0.7860*10+0 

0.7860xl0+0 

0.1000xl0+1 

0.1000xl0+1 

10'5 

0.7860x10+1 

0.7860xl0+1 

0.1000xl0+2 

0.1000xl0+2 

10“4 

0.7861xl0+2 

0.7860xl0+2 

0.1000xl0+3 

0.1000xl0+3 

Table  28 

Approximations  of  virtually 

and  nonlinear  dose 

VSD$  (t  ;it) 

safe  dose  with  k  =  4 

function 

VSDcs(t,l;*) 

Using 

Usi  ng 

Search 

Search 

Tl 

Technique  Approximation 

Technique  Approximation 

10'7 

0.4284x 10+0 

0.4284xl0+0 

0.4642xl0+0 

0.4642x10+0 

10"6 

0.9229xl0+0 

0.9229xl0+0 

0.1000xin+1 

0.1000xl0+1 

10'5 

0.1988xlOfl 

0.1988xl0+1 

0.2154xlo+1 

0.2154xi0+1 

10'4 

0.4284x 1Q+1 

0.4284*10+1 

0.4642xl0+1 

0.4642x10+1 
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Table  29 


Approximations  of  virtually  safe  dose  with  k  =  5 
and  linear  dose  function 


n 

vso, 

VSDcs(t,l;*) 

Using 
Search 
Techni que 

Approximation 

Using 

Search 

Technique 

Approximation 

10"7 

0.7634X10'1 

0.7634xl0_1 

0.1000*10+0 

0.1000xl0+0 

10-6 

0.7634*10+0 

0.7634xl0+0 

0.1000xl0+1 

0.1000xl0+1 

10-5 

0.7634xl0+1 

0.7634x10+1 

0.1000xlO+Z 

0.1000xl0+Z 

10'4 

0.7635xl0+2 

0.7634xl0+2 

0.1000xl0+3 

0.1000xl0+3 

Table  30 


Approximations  of  virtually 

safe  dose  with  k 

=  5 

and 

nonlinear  dose 

function 

VSDs 

(t;n) 

VSDcs 

(t ,l;n) 

TX 

Using 

Search 

Technique 

Approximation 

Usi  ng 
Search 
Technique 

Approximation 

10'7 

0.4242xl0+0 

0.4242xl0+0 

0.4642xl0+0 

0.4642xl0+0 

10'6 

0.9140xl0f0 

0.9140xl0+0 

0.1000xlO+1 

0.1000xl0+1 

10‘5 

0.1969xl0+1 

0.1969xl0+1 

0.2154xl0+1 

0.2154xl0+1 

10‘4 

0.4242xl0+1 

0.4242xl0+1 

0.4642xl0+1 

0.4642x10+1 

Table  31 

Approximations  of  virtually  safe  dose  with  k  =  6 
and  linear  dose  function 


V  SDs 

( t ;  it) 

VSDcs 

(t,l;*) 

Us  i  ng 

Search 

Technique 

Approximation 

Using 

Search 

Technique 

Approximation 

0.7450x10-1 

0.7450xl0_1 

0.1000xl0+0 

0.1000x10+0 

0.7450x10+0 

0.7450x10+0 

0.1000xl0+1 

0.1000x10+1 

0.7450x10+1 

0.7450x10+1 

0.1000xl0+2 

0.1000x10+2 

0.7451xl0+2 

0.7450xl0+2 

0.1000x10+3 

0.1000x10+3 

Table  3? 

Approximations  of  virtually 

safe  dose  with  k 

=  6 

and 

nonlinear  dose 

function 

VSD$ 

(t;n) 

VSD 

cs 

(t,l;n) 

Using 

Search 

Technique 

Approximation 

Usi  ng 
Search 
Technique 

Approximation 

0.4208*10+0 

0.4208xl0+0 

0.4642x10+0 

0.4642x10+0 

0.9065x10+0 

0.9065xl0+0 

0.1000x10+1 

0.1000x10+1 

0. 1953x10+1 

0. 1953x10+1 

0.2154x10+1 

0.2154xi0+1 

0.4208x10+1 

0.4208x10+1 

0.  V542x10+1 

0.4642x10+1 

The  approximations  agree  almost  always  with  the  virtually  safe  doses 


obtained  from  the  numerical  search  technique  to  at  least  four  digits. 

The  values  of  the  virtually  safe  doses  appear  to  be  changing  by  exactly  a 
factor  of  10  as  n  increases  when  the  dose  function  is  linear.  This  is 
because  it  is  increasing  by  a  factor  of  10  also.  These  approximations 
appear  to  be  much  more  accurate  than  those  for  the  mean  free  dose. 


6.2.3  Relative  Bias 


The  relative  bias  of  the  simple  virtually  safe  dose  with  respect  to 
the  cause-specific  virtually  safe  dose  is  given  by 


RB(VSDs,VSDcs)  = 


VSDs(t;*) 

VS0cs(t;n) 


The  form  of  the  specified  response's  cumulative  hazard  function  is 

the  same  as  it  is  in  the  discussion  of  the  mean  free  dose  and  is  given  by 

equation  (6.4).  The  values  of  r,  p  and  k  were  chosen  in  the  same 

-7  -6  -5 

manner.  The  relative  biases  were  calculated  for  n  =  10  ,10  ,10  and 

-4 

10  and  are  given  in  Tables  33  and  34. 

Tables  33  and  34  provide  conclusions  analogous  to  those  drawn  from 
Tables  23  and  24: 

(i)  All  relative  biases  are  negative,  indicating  that  VSDs(t;u) 
understates  VSDcs ( t ; -re ) . 

(ii)  Relative  biases  are  virtually  unaffected  by  changes  in  it.  This 
is  expected  since  u  cancels  when  the  relative  bias  is 
calculated  using  the  approximation. 


Table  33 

Relative  bias  for  virtually  safe  dose  with  linear  dose  function 


10‘7 

-0.1855 

-0.2140 

-0.2366 

-0.2550 

10  0 

-0.1855 

-0.2140 

-0.2366 

-0.2550 

10'5 

-0.1855 

-0.2140 

-0.2366 

-0.2550 

1 

O 
*— 4 

-0.1855 

-0.2140 

-0.2366 

-0.2550 

Table  34 

Relative  bias 

for  virtually 

safe  dose 

with  nonlinear 

dose  function 

10*' 

-0.06612 

-0.07712 

-0.08605 

-0.09345 

IQ"6 

-0.06612 

-0.07712 

-0.08605 

-0.09346 

10'5 

-0.06612 

-0.07712 

-0.08605 

-0.09346 

10~4 

-0.06613 

-0.07712 

-0.08605 

-0.09346 

(iii)  Relative  biases  have  greater  absolute  values  as  k  increases, 
(iv)  The  magnitude  of  the  relative  bias  is  greater  when  the  dose 
function  is  linear  than  when  it  is  nonlinear. 

The  relative  biases  for  the  virtually  safe  dose  are  about  the  same 
as  for  the  mean  free  dose.  Again  the  relative  bias  increases  with  k  and 


is  larger  when  the  dose  function  is  linear  rather  than  nonlinear.  For 
k  =  6  and  a  linear  dose  function,  the  relative  bias  of  VSD s ( t ; it)  with 
respect  to  VSDcs ( t  ,l;it)  is  -0.2550.  Thus,  VSD$(t;ii)  is  about  0.75  times 
as  large  as  VSD  (t,l;ii).  That  is,  the  virtually  safe  dose  calculated 
under  the  simple  model  is  only  0.75  times  as  large  as  the  true  virtually 
safe  dose  when  competing  risks  are  considered. 


Ui  w 


58 

7.  ESTIMATION  OF  HAZARD  FUNCTIONS  WHEN  THE  CAUSE  OF  DEATH  IS  UNCERTAIN 

The  results  in  the  preceding  section  were  developed  with  the 
assumption  that  the  specified  response  of  interest  is  observable,  e.g., 
death  with  tumor,  observable  tumor  onset,  or  death  from  tumor  (when  the 
cause  of  death  can  be  determined).  In  the  last  case,  death  from  tumor,  a 
common  problem  that  arises  is  that  the  cause  of  death  cannot  be 
determined  accurately.  This  section  describes  two  risk  assessment 
options  when  the  cause  of  death  is  unknown,  but  it  is  known  whether  an 
individual  died  with  a  tumor  or  not. 

7.1  Option  One:  All  Causes  of  Death  Combined 

Option  one  is  to  consider  the  difference  in  the  survival  function 
(or  the  corresponding  cdf)  for  a  specified  dose  d=d  versus  the  survival 
function  (or  cdf)  for  the  dose  d=0.  That  is,  compare 

Pr(T<t|d=d  )  versus  Pr(Tct|d=0)  . 

Here,  the  time  to  death  (from  all  causes  combined)  is  the  specified 

response.  If  the  zero  dose  level  is  taken  to  be  the  level  naturally 

occurring  in  the  environment,  then  Pr(T<t|d=0)  can  be  obtained  from  life 

tables  or  estimated  from  the  experimental  data.  The  probability 

Pr(T<t|d)  can  be  estimated  as  a  function  of  d  from  the  experimental 

★ 

data.  An  acceptable  dose  d=d  could  then  be  defined  as  the  dose  at  which 
Pr(T<td«d  )  exceeds  Pr(T<t  d=0)  by  no  more  than  a  specified  amount,  say 
10"6,  for  any  value  of  t  in  a  specified  time  interval.  This  approach 
assumes  that  the  dose's  effect  on  the  combination  of  all  causes  of  death 


is  of  interest.  If  only  the  dose's  effect  on  death  related  to 
carcinogenicity  is  considered  to  be  important,  then  option  two  is  more 
appropriate. 


7.2  Option  Two:  Death  Related  to  Tumor  j 

Option  two  is  derived  from  a  model  discussed  in  both  McKnight  and 
Crowley  (1984)  and  Kalbfleisch  et  al .  (1983)  but  utilizes  a  new 

definition  for  the  specified  response.  The  model  is  given  as  follows:  j 

Let 

★ 

x  -  time  to  tumor  onset  , 

x  =  time  to  death  .  I 

f 

\ 

The  hazard  rates  to  follow  use  the  notation  given  by  Table  35  as  ! 

subscripts.  j 

Table  35  ; 

Notation  for  hazard  functions  used  for  death  related  to  tumor  J 

S 

•l 

Subscript  Denotes:  j 


T 

presence  of  a  tumor 

NT 

no  tumor  present 

D 

death 

DFT 

death  from  tumor 

DCR 

death  from  competing  risk 

DRT 

death  related  to  tumor 

Define  the  hazard  rate  for  tumor  onset  to  be  , 

i 

i 


the  conditional  hazard  rate  for  death  given  that  a  tumor  has  occurred  to 
be 


X-lT(t)  -  11m  Mt<t<t+AtlT>t,T*<t)  , 
°|T  At+O  « 


and  the  conditional  hazard  rate  for  death  given  that  no  tumor  has 
occurred  to  be 


xnl„T(t)  =  11m  Pr(to.<t+At|T>t,x  >t)  _ 

U|  Nl  ...  At 


At+0 


Also  let 


J  =  < 


1  if  the  death  is  due  to  tumor 

^2  if  the  death  is  due  to  competing  risks. 

Then  the  hazard  rate  Xp|T(t)  can  be  partitioned  as  follows: 

x  ,  (t)  *  11m  Pr(t<T<tmAt,J-l|T>t,m*<t) 
t  At+o  ^ 


+  lim  Pr(t<T<t+At,Js2  T>t,t  <t) 

At+0  At 

XDFT J ^  +  XDCR J 1  ^  ’  (7 

where  the  last  two  hazard  rates  are  the  hazards  corresponding  to  death 
from  tumor  given  that  a  tumor  has  occurred  and  death  from  a  competing 
risk  given  that  a  tumor  has  occurred,  respectively.  Notice  that  the 
conditional  hazard  rates  do  not  depend  on  the  time  of  the  tumor  onset 
rather  they  depend  only  on  the  existence  of  the  tumor  prior  to  time  t. 
Mcknight  and  Crowley  (1984)  discuss  this  situation  with  respect  to 
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Xp|y(t).  If  this  hazard  rate  is  written  as  being  conditional  on  the 
time  of  tumor  onset,  then  it  is  not  identifiable.  Therefore,  X^^t)  is 

regarded  as  a  marginal  hazard  function  that  averages  the  risk  over  the 

★ 

distribution  of  x  . 

If  the  cause  of  death  can  be  determined,  then  \j(t),  XDCR|-]-(t), 
^DFTjT  anc*  xd|nT^  can  ke  sstimated  provided  interim  sacrifices  are 
made  (see  Mcknight  and  Crowley  (1984)).  Kalbfleisch  et  al .  (1983)  give 
the  likelihoods  for  this  case.  However,  if  the  cause  of  death  cannot  be 

T(t)  are  not 

identifiable  and  only  j ^.( t )  can  be  estimated. 

Assuming  that  the  cause  of  death  cannot  be  determined,  the  focus  is 
now  shifted  to  a  new  specified  response  -  death  related  to  tumor.  This 
response  has  meaning  only  for  the  case  where  the  hazard  rate  in  the 
presence  of  a  tumor  exceeds  the  hazard  rate  in  the  absence  of  a  tumor, 
death  related  to  a  tumor  is  defined  to  be  a  death  caused  by  the  increase 
in  Xp  |  y(  t)  relative  to  X0  j  T: ) .  This  implies  that  the  hazard  rate  for 
the  specified  response  is 

xdrt|t^  =  xd|t^  "  xd|nt^  *  (7*2^ 

A  death  related  to  tumor  may  be  a  death  caused  by  a  tumor  directly  or  a 
death  resulting  from  an  increase  in  the  hazard  rate  for  a  competing 
risk.  Since  a  death  related  to  tumor  is  not  necessarily  clinically 
differentiable  from  other  deaths  with  tumor,  it  is  important  that  the 

^(t)  does  not  require  that  such 


estimation  procedure  for  XpR^ 


determined,  then  the  hazard  functions  and  X^ 


differentiations  be  made. 
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If  is  replaced  by  its  equivalent  expression  in  (7.1),  then 

xdrt|t^  =  xdft|t^  +  xdcr|t^  "  xd|nt^  • 

which  futher  emphasizes  that  a  death  related  to  tumor  is  not  just  a  death 
from  tumor  but  refers  to  the  increase  in  all  causes  of  death  in  the 
presence  of  a  tumor  relative  to  all  causes  of  death  in  the  absence  of  a 
tumor. 

To  estimate  the  desired  risk  characterizations,  it  is  necessary  to 
first  estimate  Xy(t),  XD|y(t)  and  Ny(t).  Tbis  *s  done  using  the 
likelihoods  from  Kalbfleisch  et  al.  1983)  which  simplify  when  the  cause 
of  death  is  unknown.  The  resulting  likelihood  contributions  are 

LI  =  XD|NT(t)  exp[‘^0  (XD|NT(u)  +  Vu))  du]  • 

Ln  *  exp[-/b  (^0|NT(u)  +  Xj(u))  du]  , 

L 1 1 1  =  Jg  XT(u)*exp[-/J  (XT(v)  +  XD | NT( v) )  dv  -  Jb  XD|T( v)  dv]du, 
and 

LIV  =  /o  VU)Xexp^0  +  XD|NT(V)^  dv^ 

x  XD|T^t^XeXpl-"^u  XD|T^  dv^  du  * 

where  Table  36  gives  the  situation  governing  each  of  the  contributions. 
Using  these  contributions,  the  maximum  likelihood  estimates  \y(t), 

\D|y(t)  and  X 


0|Ny(t)  can  be  obtained  and  XDRyi 


y ( t )  is  obtained  by 
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Table  36 

Situations  for  the  likelihood  contributions  when 
the  cause  of  death  is  uncertain 


Contribution 


Situation 


death  at  time  t  with  no  tumor  present 
sacrifice  at  time  t  with  no  tumor  present 
sacrifice  at  time  t  with  a  tumor  present 
death  at  time  t  with  a  tumor  present 


applying  equation  (7.2).  Thus,  the  time-to-response  probabilities  and 
other  risk  characterizations  for  the  time  to  death  related  to  tumor  can 
be  estimated. 

The  dependence  of  the  hazard  functions  on  dose  has  been  suppressed 
in  this  discussion.  A  natural  experimental  situation  would  be  to 
estimate  the  probability  of  death  related  to  tumor  at  dose  d  and  compare 
this  to  the  probability  of  death  related  to  tumor  at  dose  0.  Also,  other 
risk  characterizations  based  on  dose,  such  as  mean  free  dose  and 
virtually  safe  dose,  can  be  estimated. 

These  two  options  provide  the  researcher  and  risk  assessor  with 
alternatives  when  the  cause  of  death  is  unknown.  Both  options  require 
that  a  new  specified  response  be  determined.  The  first  option  gives  risk 
characteri zations  when  death  from  all  causes  combined  is  the  specified 
response.  The  second  option  gives  a  method  for  estimating  risk 
characteri zati ons  when  death  related  to  tumor  is  the  specified  response. 


•  .vV* 
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8.  CONCLUSIONS 

This  dissertation  has  demonstrated  a  method  for  improving  risk 
characterizations  from  a  simple  model  by  incorporating  competing  risks 
via  the  cause-specific  model.  Two  aspects  of  this  method  arise  which 
make  the  incorporation  of  the  cause-specific  model  desirable;  they  are 
increased  accuracy  and  ease  of  implementation. 

The  cause-specific  model  give  mure  accurate  representations  of  risk 
characteristics  than  the  simple  model  in  two  ways. 

(i)  The  cause-specific  model  reflects  the  possibility  of  a 

competing  response  occurring  before  the  specified  response 
whereas  the  simple  model  does  not. 

(ii)  The  cause-specific  model  implies  more  relevant  risk 
character!'  zations. 

The  difference  between  the  two  models  was  discussed  in  terms  of 
relative  bias.  In  Sections  5  and  6,  bounds  are  obtained  for  various  risk 
characterizations  and  examples  are  given  that  demonstrate  how  the  simple 
model  overstates  the  effect  of  the  carcinogen.  If  the  risk 
characterization  is  a  time-to-response  probability,  the  simple  model  can 
overstate  the  probability  by  as  much  as  800%  depending  on  latency  period, 
time,  and  parameter  values  of  the  specified  response's  hazard  function. 
However,  most  overstatements  are  of  the  order  of  50%  -  60%  for  middle 
values  of  time,  latency  period  and  parameter  values. 

Additionally,  if  the  data  are  generated  from  a  cause-specific  model 
with  a  latency  period  but  fit  by  a  simple  model  without  a  latency  period, 
then  the  overstatement  of  the  carcinogenic  effect  is  usually  increased. 
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This  was  demonstrated  by  a  simulation  study  in  Section  5.  It  shows  that 
the  simple  time-to-response  probability  may  be  as  much  as  1400%  larger 
than  the  true  time-to-response  probability  under  the  cause-specific 
model . 

The  magnitude  of  the  relative  bias  for  mean  response  free  period  was 
not  as  large  as  it  was  in  the  case  of  time-to-response  probability.  In 
fact,  if  the  probability  of  the  specified  response  occurring  before  time 
t  is  of  the-  order  10“®  then  the  relative  bias  for  mean  response  free 
period  is  almost  negligible;  however,  if  the  risk  characterization  is 
given  in  terms  of  dose,  the  relative  bias  may  be  larger.  For  example,  if 
the  dose  function  is  linear,  the  virtually  safe  dose  under  the  simple 
mmodel  may  only  be  75%  of  the  virtually  safe  dose  under  the 
cause-specific  model.  The  same  order  of  magnitude  holds  for  the  relative 
bias  of  the  mean  free  dose.  This  percentage  increases  to  about  90%  if 
the  dose  function  is  nonlinear  (e.g.,  cubic). 

The  risk  characteristics  under  consideration  have  different 
representations  under  the  simple  and  cause-specific  models.  By  using  the 
new  definitions  of  the  risk  characterization  under  the  cause-specific 
model  the  researcher  is  required  to  be  more  precise  in  his  definition  of 
the  specified  response  and  the  health  effects  to  be  measured.  These  risk 
characteristics  under  the  cause-specific  model  are  easily  interpreted. 

Accuracy  is  just  one  advantage  of  the  cause-specific  model; 
implementation  is  another.  As  shown  in  Section  3,  the  relationship 
between  the  cause-specific  and  simple  likelihood  functions  allows  the  use 
of  existing  simple  model  maximum  likelihood  estimates  to  develop 
cause-specific  maximum  likelihood  estimates.  Also,  the  approximation  for 


the  competing  risks'  hazard  function,  X2(t),  given  in  Section  4  is  both 
reasonable  and  convenient.  Thus,  the  procedure  given  in  Sections  3  and  4 
for  incorporating  the  cause-specific  model  can  be  implemented  using 
either  the  experimental  data  alone  or  the  simple  model  estimates  and  a 
life  table  approximation  for  the  competing  risks'  hazard  function. 

Finally,  when  the  cause  of  death  is  uncertain,  the  framework  of  the 
cause-specific  model  allows  the  estimation  of  time-to-response 
probabilities  for  the  specified  response,  death  related  to  tumor.  Using 
this  new  response  is  reasonable  and  may  be  more  relevant  to  the  risk 
assessment  than  defining  the  response  to  be  death  from  tumor  when  the 
cause  of  death  is  unknown  or  subjectively  assigned. 

In  the  area  of  quantitative  cancer  risk  assessment  many  modeling 
problems  occur  such  as  low  dose  extrapolation,  interspecies  extrapolation 
and  the  incorporation  of  competing  risks.  These  are,  by  no  means,  easy 
problems  to  solve.  They  should  each  be  considered  by  risk  managers,  risk 
assessors  and  researchers  with  the  attitude  that  accuracy  in  all  aspects 
of  the  modeling  process  is  important.  This  dissertation  provides  these 
people  with  the  necessary  information  to  decide  if  the  simple  model  is 
sufficient  in  a  particular  situation,  and  if  it  is  not,  it  provides  a 
method  for  improving  the  risk  characterizations  by  utilizing  existing 
results  without  setting  up  new  experiments. 
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